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Introduction: 

STARS is a multicenter NASA project to determine the feasibility of using space-based 
assets, such as the Tracking and Data Relay Satellite System (TDRSS) and Global 
Positioning System (GPS), to increase flexibility (e.g. increase the number of possible 
launch locations and manage simultaneous operations) and to reduce operational costs by 
decreasing the need for ground-based range assets and infrastructure. The STARS 
project includes two major systems: the Range Safety and Range User systems. The 
latter system uses broadband communications (125 kbps to 500 kbps) for voice, video, 
and vehicle/payload data. Flight Demonstration #1 revealed the need to increase the data 
rate of the Range User system. During Flight Demo #2, a Ku-band antenna will generate 
a higher data rate and will be designed with an embedded pointing algorithm to guarantee 
that the antenna is pointed directly at TDRS. This algorithm will utilize the onboard 
position and attitude data to point the antenna to TDRS within a 2-degree full-angle 
beamwidth. This report investigates how errors in aircraft position and attitude, along 
with errors in satellite position, propagate into the overall pointing vector. 

Explanation of Variables between the Aircraft and TDRS: 

The antenna on the aircraft must have a direct line-of-sight to TDRS, which is 
represented by the green arrow in figure 1 . The vector is defined by azimuth (a) and 
elevation (e) angles in the antenna platform coordinate system. The antenna’s actual 
pointing vector, represented by the blue arrow and a' and s', might differ from the desired 
vector pointing to TDRS. The 8 angle is the total difference between the TDRS and 
antenna directions. This angle must have less than a 2 degree full-angle beamwidth in 
order for the aircraft’s antenna to view TDRS. 
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Figure 1. Antenna and TDRS Pointing Vectors 


Coordinate Transformations: 

TDRSS is nearly geostationary, in a low-inclination orbit with known position 
represented as latitude, longitude, and altitude. The aircraft’s position is represented in 
the Earth-Centered Earth-fixed (ECEF) coordinate system and its attitude represented as 
heading, pitch, and roll. The ECEF coordinate system has three axes labeled e, f, and g. 
The g-axis points from the center of the earth to the north pole (along the axis of 
rotation); the e-axis points from the center of the earth to the equator at the prime 
meridian; and the f-axis is derived from the e- and g-axes to be orthogonal, pointing into 
the Indian Ocean at the equator. 


The Dennis.c method is used to convert the TDRSS position from latititude, longitude, 
and altitude into ECEF coordinates. This method first defines r e / g as the radius vector to 
TDRS in the ECEF system. The latitude and longitude of TDRS (<pr and At) are known, 
along with the distance from TDRS to Earth center for geosynchronous range ( Rt ), which 
can be calculated using Kepler’s Third Law. 
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These three values represent the E, F, and G coordinates of TDRSS, respectively. The 
aircraft’s E, F, and G coordinates were included in the flight data for flight 6. The 
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aircraft’s ECEF position vector is subtracted from the TDRS position vector to form the 
first pointing vector P. 

This vector must now be transformed from the ECEF coordinate system to the Local 
Horizontal Coordinate System. This coordinate system is also called the earth-fixed axis 
or the geographic coordinate system. One axis points to the north, one to the east, and 
one axis is either up or down. For this application we will be using the north-east-up 
method (NEU). The local horizontal system is related to the ECEF system through 
latitude, longitude, and altitude. The specific method for the NEU transformation is , 
provided below, and results in the second pointing vector, P’: 
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Once the pointing vector is in the Local Horizontal Coordinate System, the corresponding 
matrix must be rotated about the NEU coordinates in order to take into account the 
attitude of the aircraft. After the rotation, the new coordinates are in the Body- Axis 
coordinate system, also known as the body platform coordinate system. It is a right- 
handed system with x-axis forward, y-axis to the right, and z-axis pointed down. This 
coordinate system is related to the local horizontal system through the heading, pitch, and 
roll angles of the aircraft. 


The specific method for the transformation is included below, and results in the third 
pointing vector, P”. The order of rotation is of particular importance: 
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The final azimuth and elevation angles can be derived from the X, Y, and Z components 
of the P” vector. 
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Azimuth : a = tan 


-i 


r Y 

p" 

\ r x J 


Elevation : s = tan 
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The following section will explain how these values were found to be accurate. 

Comparing the Azimuth and Elevation Results with Phi and Theta Aircraft Data: 

To ensure the validity of the series of transformations and rotations performed in the 
previous section, the equations were checked using actual data from Flight 
Demonstration #1 . The data included pointing information in another coordinate system 
from CLASS, and this data was compared to the data derived in the previous section. 
The phi and theta data used in the system from CLASS are defined identically to the 
STARS Flight Demo #1 Report. Specifically: 

• theta = 0°, phi = any value is the nose of the F-l 5 

• theta = 1 80°, phi = any value is the tail of the F-l 5 

• theta = 90°, phi = 0° is straight up 

• theta = 90°, phi = 90° is out the starboard wing 

• theta = 90°, phi = 1 80° is straight down 

• theta = 90°, phi = 270° is out the port wing 

The NED coordinate system was used initially to define phi and theta. The phi and theta 
coordinates from the flight data were used to convert to azimuth and elevation using the 
following: 



Additionally, 

x = R'cos(0) _y = z'sm(^) = Rsm(#)sin(^) 

z'=7?'sin(^) z = z'cos(^) = i?'sin(#)cos(#>) 
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Elevation(deg) Azimuth(deg) 


The final azimuth and elevation calculations can then be derived: 
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The two separate calculations for azimuth and elevation were then plotted against each 
other: 


Algorithm Azimuth vs. Phi/Theta Derived Azimuth for Flight 6 



Algorithm Elevation vs. Phi/Theta Derived Azimuth for Flight 6 
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Algorithm Azimuth vs. Phi/Theta Derived Azimuth for Flight 7 



Algorithm Elevation vs. Phi/Theta Derived Azimuth for Flight 7 



Errors Added into the Algorithm: 

After confirming the azimuth and elevation values, errors were added into the position 
and attitude of the aircraft at each time increment. The “perturbed” values of position 
and attitude were subtracted from the original values to form a new “perturbed” value of 
azimuth and elevation. In addition, the overall angle between the antenna’s pointing 
vector and the TDRS pointing vector, called delta, was calculated using the spherical law 
of cosines, or: 

TDRS ant ~ (l)(l)cos£ 

= sin e sin s' + cos e cos s' cos(a - a') 

The variables in this equation have been defined previously. Since all variables have 
known values, the delta angle can be found simply using trigonometry. The azimuth and 
elevation delta values were then plotted against time to reveal the effects of the errors. 

Results: 

Aircraft position errors seemed to have a much reduced effect on the data than did 
changes in aircraft attitude or TDRS position. A 1 -degree change in heading, pitch, or 
roll created significant differences in azimuth and elevation delta values, but even a 
10,000 meter change in position had little effect. Change in TDRS position had an even 
greater effect, with a 0.05 degree change in latitude or longitude pushing the overall delta 
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value beyond the 2 degree full-angle beamwidth (see plot below). The numbers next to 
the title of the plot signify the error of terms in the following order: x-position / y- 
position / z-position / heading-angle / pitch-angle / roll-angle. Also included are the 
minimum and maximum delta and the times at which they occurred, the mean delta, and 
the standard deviation of the delta. The plot below includes a dotted line at 2-degrees to 
distinguish between valid and invalid data. It was taken from flight 6 data; flight 7 data 
was taken at a higher frequency and had a much greater change in the magnitude of the 
delta values. 


[Delta Beamwidth Angle: 0/0/0/0/5/0 
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Summary of Results: 
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E (aircraft) 

-10km — ► 10km 

-0.1 0.1 

-0.1 ^0.1 

0.1 

6 

miH 

F (aircraft) 

-10km —> 10km 


^ | 

0.06 

6 

lOUHi 

G (aircraft) 
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14 

H (heading) 

-20° -> 20° 

-10 ^ 10 

-28 -*• 28 

19.2 ±0.1 

6 
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-12 — >■ 12 


6 
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-17 —*■ 17 

BEE 

6 
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TDRS Lat. 

-2° -> 2° 

-2 — ► 2 

-2.5 2.5 

2.1 

6 

8 

TDRS Long. 

-2° -> 2° 

-2 — >■ 2 

-2.5 -»• 2.5 

2.1 

7 

10 

E (aircraft) 

-10km — ► 10km 

-0.1 0.1 

-8.9 7.6 

0.1 

7 
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F (aircraft) 

-10km — > 10km 

-0.06 -»• 0.06 

-3.9 3.4 

0.06 

7 

10 

G (aircraft) 

-10km — *• 10km 

-0.05 -+ 0.05 

-3.0 3.3 

0.05 

7 

14 

H (heading) 

-20° 20° 

-19.2 — ► 19.2 

-175 ^ 180 

19.2 ± 0.1 

7 

14 

P (pitch) 

-20° -»• 20° 

-20 -+ 20 

-171 -> 169 

5.33 — »■ 20 
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R (roll) 

-20° -*• 20° 

-20 20 
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00 

0 

1 

00 

o 

0.01 -*■ 20 
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TDRS Lat. 

; -i°-»i° 

-2.1 — ► 2.1 

most within ± 10°, 
some up to ± 50° 

2.1 

7 

2 

TDRS Long. 

-i°-> r 

-2.1 — *• 2.1 

most within ± 1 0°, 
some up to ± 50° 
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